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Abstract. Experimental data of arrays "Andyrchy" 
and "Carpet-2" of Baksan Neutrino Observatory (In- 
stitute for Nuclear Research), obtained in the regime 
of a detection of the single cosmic-ray component, are 
used for a search of the bursts of cosmic gamma rays 
from evaporating primordial black holes. Different 
theoretical models of the evaporation process are 
used for the analysis. Distributions of the counting 
rate fluctuations on both arrays agree with the 
expectations from the cosmic ray background. The 
new constraints on the concentration of evaporating 
primordial black holes in the local region of Galaxy 
are obtained. The comparison of the results of 
different experiments is given. 
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I. Introduction 

It has been argued recently |1| that the photon flux 
calculation from an evaporating black hole (BH) given in 
IJl is the most reliable, and the photospheric and chro- 
mospheric effects considered in O, ||4| are negligibly 
small. It is clear, however (and we try to show this in 
this section), that scenarios of BH evolution, in which 
interactions between emitted particles (and even some 
kind of the thermal atmosphere around the BH) exist, 
cannot, in general, be discarded. 

In Sees, [n] and Ull] we study the detection of BH 
evaporation using three models (considering them on 
an equal footing). We do not insist on the validity of 
concrete chromospheric models; our aim is to demon- 
strate the sensitivity of the experimental method used in 
the present paper (which is based on the registration of 
evaporated photons with energies ~ 10 GeV) to a form 
of the BH photon spectrum, having in mind that its right 
shape may be different from the "canonical" one. 

The famous Hawking 's result ||5l according to which 
a BH will emit black body radiation, precisely corre- 
sponding to the temperature Th — /^ornsH (^^^^"^ '^R 
is the dimensionless Rindler temperature, Tj^ = t-) had 
been obtained in a semiclassical approximation. It had 



been assumed, in particular, that back-reaction effects 
are effectively small and do not influence the classical 
collapse geometry. It is well known that Hawking 's 
derivation of the evaporation spectrum is based on 
information loss phenomenon |6| and, correspondingly, 
on the possibility of an nonunitary evolution of the BH. 
This is a direct consequence of the strict locality of 
quantum field theory which leads to an independence 
of the external and internal Hilbert states (i.e., states 
located outside and inside horizon, respectively) and, as 
a result, to an appearing of the density matrix for the 
BH decay rather than a pure state. 

It had been argued in many works, however (see, 
e.g., |121> EIj ID), that the semiclassical approximation 
used in f5| is too crude: generically, one cannot regard 
outgoing particle as a Hawking particle that was not af- 
fected by the ingoing matter Moreover, the gravitational 
interactions of the Hawking radiation (when it already 
separates from the horizon) with infalling particles can 
be strong enough for a deforming the metric and for a 
breakdown of locality. It is known that the dynamics of 
strong gravitational physics is inherently nonlocal, and 
it is natural to assume that just this nonlocality may 
be the resolution of the above-mentioned paradox of 
information loss |9|. 

An approximate (effective) accounting of the back- 
reaction effects is realized in the approach in which 
the boundary condition for all fields on a surface a 
few Planck distances away from horizon is postulated. 
The examples are the "brick wall" model |10|, the 
"stretched horizon" model ifTTl and the "bounce" model 
lfT2l . These models are based on an assumption that 
strong gravitational interactions among the field quanta 
effectively form this barrier between the horizon and the 
wall. 

It is argued in these works that the process of forma- 
tion and evaporation of a BH, as viewed by a distant 
observer, can be described entirely within the context 
of standard quantum theory, with unitary S'-matrix and 
pure quantum states. 

An idea of the stretched horizon (SH) was proposed 
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many years ago ifTSll but at that times it was considered 
as a useful mathematical construction rather than a phys- 
ical object equipped with some microphysical degrees of 
freedom. One should stress, however, that the approach 
using a concept of the SH is appropriate only for an 
external observer, not for a free falling one. The main 
point of the approach is that the SH can, through its 
degrees of freedom, absorb, thermalize and emit any 
quantum mechanical information falling into the BH, 
without any information loss ifTTI . 

The local proper temperature at the SH, T,, is related 
to the temperature measured by distant observers (which 
is, by definition, the Hawking temperature T/f), by the 
connection 



dt . 
d^ 



(1) 



where dr/dt is the time dilation factor, connecting the 
time intervals at the SH with the intervals of coordinate 
time t, which is given by 

rfr _ ph ^ Ph_ 
dt ~ UIG ^ Th ' 



(2) 



AMG 

Here, ph is the distance of SH from the horizon, Vh is 
the gravitational radius. In the "standard" case, when 
Ph ~ GMp = Ip {Mp is the Planck mass), one has 

^ ^ ^ (3) 
dt M ' ^ ' 

and, if Tg ^ Mp (independently of the size or mass of 

the BH), one obtains 



Th 
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It is very important that the SH is assumed to be 
in a thermal equilibrium with the surrounding material 
during most of the evaporation. A distant observer would 
estimate the number of particles emitted per unit time 
which is proportional to a product of the BH area and 
the time dilation factor ITTI . 
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(if all these particles go to infinity). On the other hand, 
the number per unit time of particles that actually 
emerged to infinity is 

dN L 1 
~dt ^ 'e< 



■typ 



M 



(6) 



(i?typ is the typical energy of the emitted particles, 
Etyp Th ^ Mf>/M, L is the BH luminosity, 
L ~ 1/Af2). It follows from © and © that most of 
the particles emitted from the SH do not go to infinity. 
This gives rise to a thermal atmosphere above the SH 
(due to repeated interactions of emitted particles with 
the SH and with each other). 

One should stress that, in these models, BH decays 
as a pure quantum mechanical state, so, a spectrum of 
the emitted radiation does not need to have a precise 
thermal form predicted in ||5l, at least at final stages of 
the evaporation, when BH mass is small. 



In general case, the temperature of the SH, Ts, as well 
as the distance ph, depend on microphysics. Correspond- 
ingly, the maximum value of the Hawking temperature 
Th, and a form of the spectrum of evaporated particles 
are model dependent. Evidently, this model dependence 
is especially large at final stages of the evaporation when 
a radius of the BH, r^, is comparable with ph- In this 
case the centrifugal barrier becomes inefficient and the 
thermal atmosphere around the SH gradually becomes 
open for a space surrounding the BH. The experimental 
signature may be similar with those predicted in chro- 
mospheric models if the maximum Hawking temperature 
Th is much smaller than Mp, as is predicted in some 
scenarios. For example, such a situation is possible in 
string models. In these models, the SH is placed at a 
distance ^ Ig from the event horizon {Is is the string 
scale). The local proper temperature at the SH is equal to 
|fT4l Ts — 1/2tiIs. We suppose, that the string coupling, 
~ Gls, is extremely small, so that the Planck and 
string scales are well separated. When, in a course of the 
evaporation, the Hawking temperature Th approaches 
Ts (which is the Hagedorn temperature of string theory), 
the radius of the BH becomes equal to Ig- In this point 
the BH can transform into a higher-entropy string state. 
The possibility of such transformation was discussed in 
many works (see [151 and references therein). 

It is essential also that the ideas of the string-BH 
correspondence |16| and of identifying the states of a 
BH with the highly excited states of a fundamental string 
( lfT4l . ifTTl ) help us to understand the physical nature 
of the BH entropy (by other words, the nature of the 
"internal states of a BH", or, "the degrees of freedom of 
the stretched horizon"). 

It is important to note that strings are ideally suited 
for interactions with the SH due to a spread in the 
transverse space (during approaching the SH) and due to 
Regge behavior of string scattering amplitudes and their 
growing cross sections |18|. 

The most important conclusion from this scenario is 
that the resulting, weakly coupled, string decays at the 
Hagedorn temperature, and the momenta of the emitted 
particles never exceed Ts (which can be as small as 0.1 
TeV). 

II. The experiment 

Up to now, the search for the bursts of high-energy 
gamma rays which are generated at the final stages of 
PBH evaporation was carried out in ground-based exper- 
iments detecting the EASs on several shower arrays llT9l . 
ESI, lEl, EH, ED and on the Whipple Cherenkov tele- 
scope ED- Because of high energy detection threshold, 
the interpretation of the results of these experiments can 
be performed only within the framework of the evap- 
oration model without a chromosphere El (hereafter, 
MW90). The duration of the high-energy burst predicted 
by chromospheric evaporation models is too short, much 
shorter than the detection dead time in these experiments 
(the duration of burst tb is, by definition, the time interval 
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Fig. 1: Time-integrated photon spectra from a BH with 
initial Hawking temperature Th = 10 TeV. Such a BH 
has about 0.5 s until the full evaporation. 



during which 99% of photons that can be detected by 
the array evaporate). 

For the direct search of PBH gamma ray bursts within 
the framework of the chromospheric models of BH 
evaporation, another technique has to be applied, which 
is sensitive to photons with rather small energies, namely 
of order of few GeV. Shower particles generated in 
atmosphere by photons with energy ~ 10 — 100 GeV 
are strongly absorbed before reaching the detector level, 
so, the average number of signals in the detector module 
is smaller than 1. Therefore, in this energy range PBH 
bursts can be sought for by operating with the modules 
in single particle mode, that is, by measuring the single 
particle counting rate of the individual modules. The 
primary arrival directions of photons are not measured, 
and bursts can be detected only as spikes (short-time 
increases) of the cosmic ray counting rate. The effective 
energy of the primary gamma-ray photons detected by 
this method depends mainly on the altitude of the array 
above sea level. This technique was employed earlier to 
seek for cosmic gamma-ray bursts with energies exceed- 
ing several GeV 125], ||26l, fTT\, f28\. First constraints 
on the number density of evaporating PBHs within 
chomospheric framework were obtained in work ||29l 
for the evaporation models by Heckler f3\ (hereafter, 
H97) and Daghigh and Kapusta [4] (hereafter, DK02). 
It should be noted that for the PBH search within the 
framework of model MW90, the same technique still can 
be applied. The comparison of three evaporation models 
considered can be found in ll30l . and time-integrated 
photon spectra for particular value of Hawking temper- 
ature are shown in Fig. [T] 

The present experiment was performed on Andyrchy 
and Carpet-2 arrays of Baksan neutrino observatory of 



INR RAS. Both arrays are situated close to each other 
(horizontal distance between them is about 1 km), but 
their altitudes above sea level are different: Andyrchy 
- 2060 m (atmospheric depth - 800 g cm~^). Carpet - 
1700 m (atmospheric depth - 840 g cm~^). 

Andyrchy consists of 37 scintillation detectors. The 
area of each detector is 1 m-^. To detect a single cosmic- 
ray component, the total counting rate of all of the de- 
tectors is measured each second. The search for gamma- 
ray bursts using this technique is carried out against the 
high cosmic -ray background (uj ~ 11440 s~^), which 
requires the highly stable and reliable performance of all 
equipment. The monitoring is realized through simulta- 
neous measurements (with 1-s acquisition rate) of the 
counting rates in the four parts of the array comprising 
10, 9, 9, and 9 detectors, respectively. An use of such 
information allows to exclude 1 -second intervals with 
unreasonably large deviations of counting rate between 
the parts of the array which could come from faultiness 
of individual detectors or registration channels, faultiness 
of counting rate measurement channels of separate parts 
of the array or non-synchronous impulsive electromag- 
netic disturbance in detectors or signal cables. A detailed 
description of the array and its operating parameters is 
given in ||3TI . 

Carpet-2 array ll32ll . ll33l consists of the actual Carpet 
facility (400 liquid scintillation detectors covering area 
of 196 m^), six remote points (RP; 108 detectors of the 
same kind with total area 54 m^), and the muon detector 
(175 detectors of the type used in Andyrchy array placed 
in the underground tunnel). 

The detection probabilities P{E^. 0) of the secondary 
particles produced by the primary photons with en- 
ergy falling on the arrays at zenith angle 6 were 
determined by simulating electromagnetic cascades in 
the atmosphere and the detectors. The effective energy 
of gamma rays registered by Andyrchy and RP is 8 
GeV, for Carpet and muon detector this energy is, 
correspondingly, 200 GeV and 2 TeV. Because of this, 
only data from Andyrchy and RP was actually used in 
this work. 

The total number of gamma rays which can be de- 
tected by the array is given by 



oo 

j dE^P{E^,e)dN^/dE^. (!) 



It depends on the time ti left until the end of PBH 
evaporation. Here, dNj/dE-y is the spectrum of photons 
emitted by the PBH during the same time interval ti. 
Let a PBH be located at distance r from the array with 
area S and be seen from it at zenith angle 9. The mean 
number of gamma-rays detected by the array is then 



n{9, r) 



N^{0,ti)S{0) 

47^7-2 



(8) 



which will give an excess over the average counting rate 
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Fig. 2: Full number of gamma particles that can be 
detected by the Andyrchy array for 6* = 0°, as a function 
of time until the end of PBH evaporation, for different 
evaporation models: 1 - DK02, 2 - H97, 3 - MW90. 

u) in units of root mean square deviation: 

The time interval At which should be used to search 
for bursts of gamma rays from evaporating PBHs in 
the particular experiment depends both on the burst 
duration and on its time profile. Burst duration for 
Andyrchy and RP {Eth ^ 8 GeV) exceeds 10"' s for all 
considered evaporation models. It should be noted that 
the search of long bursts (tens of seconds and more) 
in ground-based experiment, especially when working 
in single-particle mode, is difficult because of short- 
time variations of cosmic ray intensity [26], 128|. Fig. 
|2] shows, for considered models, the dependencies of 
total number of gamma rays that can be detected by 
Andyrchy array on the time ti. It is seen that even for a 
rather strong increase of ti (from 1 s to 20 s) the number 
of gamma particles increases only by a factor of 2 to 4 
(depending on the model). At the same time, the number 
of background events increases proportional to the time 
interval. The consequence of this consideration is that in 
our experiment, the search is most effective for 1 -second 
time intervals. 

Deviations in the array counting rate lasting At < 1 
s are sought using the parameter Fi that is equal to the 
deviation (measured in units of the Poisson sigma) of 
the number of counts ki during the i-th second of a 15- 
min interval from the average number of counts during 
this interval: Fi — {ki — k)/Vk. Since variations in the 
cosmic -ray intensity over a time of 15 min are negligible 
in the first approximation and the average counting rate 
is fairly high, one can expect that the parameter Fi has 
a Gaussian distribution with the zero mean value V = 
and unit standard deviation a = 1.0. 

The actual experimental data mostly confirms this 
assumption. For the RP data, no excesses were found 
larger than 7 standard deviations. For Andyrchy, the only 
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Fig. 3: Maximum distance from which the evaporating 
PBH can be detected. Solid lines are for Andyrchy 
(number of particles in this case is n = 640), dotted 
lines are for RP (n = 918); models used are: 1 - DK02, 
2 - H97, 3 - MW90. 



event with a large (7.9a) deviation was detected on April 
17, 2002 at 17:31:29 UT , all other data is well fitted 
by a Gaussian distribution with maximum deviations of 
6(7. If we assume that the 7.9(t Andyrchy event is caused 
by evaporating PBH, then, taking into account that RP 
is examining the same region in the sky, the excess in 
RP counting rate should be in the region (6.6 — 7.7)(t 
for chromospheric models, or in the region (8.2 — 8.5)(t 
for MW90 model. However, for this 15 -min interval, no 
excess larger than 3.2(7 was detected by RP, so this single 
Andyrchy event cannot be explained by the evaporating 
PBH. 

One can estimate the maximum distance from which 
a PBH can be seen using the formula [see Fig. [3] 

ao) 

III. Results and conclusions 

Generally, the effective volume surveyed by the array 
is calculated using the formula 

oo 

V^S = J dnj drr^F{n,n{e,r)). (11) 



Here, F{n, n{d, r)) is the Poisson probability to register 
n or more events having the average n{9, r) determined 
by Eq. Taking into account that no excesses above 
6(7 were detected on Andyrchy, we put n = 6\/q = 640, 
and the effective volume for this array is, depending on 
the model: 1.2 xlO^^" pc^ 5.6x10"^° pc^, 8.1x10"" 
pc^ for H97, DK02 and MW90, correspondingly. For 
RP, no excesses above 7a are discovered (in this case 
u} = 17200, so, n = 7Vq = 918), and the effective 
volume is: 7.6 x 10"" pc^, 4.1 x lO^^*' pc^, 7.1 x 10"" 
pc^ for H97, DK02 and MW90, correspondingly. 
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The number of bursts detected over the total observa- 
tion time T can be represented as 

N - PpbhTFeff , (12) 

where ppbh is the number density of evaporating PBHs. 
Assuming that evaporating PBHs are distributed uni- 
formly in the local region of the Galaxy and taking into 
account that both arrays survey the same sky region at 
different time, one can calculate the upper limit piim 
on the number density of evaporating PBHs at the 99% 
confidence level using the formula 

Pnr,,^A.6/{VATA + VRpTRp), (13) 

with full observation time Ta = 6.27 yr for Andyrchy 
and Trp = 2.34 yr for RP. Substituting the effective 
volumes with values calculated above, we finally obtain 
for piim: 10^ pc-3 yr-i, 5 x 10^ pc'^ yr"!, 6.8 x 10^ 
pc~^ yr^^ for evaporation models DK02, H97 and 
MW90, respectively. 

One can note, that in a case of the non-chromospheric 
model MW90, the limit obtained in this work is worse 
than the limit from the Andyrchy array obtained using 
method of EAS detection |23|. However because the 
effective energies of detected photons differ by several 
orders of magnitude, these Umits should be regarded as 
complementary. For the case of DK02 and H97 models, 
the limit has been significantly improved compared to 
our previous work ||29]| . 
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